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Abstract 

In the paper the dynamic localization of charged particle (electron) in a quan- 
tum wire under the external non-uniform time-dependent electric field is consid- 
ered. The electrons are trapped in a deep 'dynamic' quantum wells which are 
the result of specific features of the potential imposed on 2D electron gas: the 
scale of spatial nonuniformity is much smaller then the electron mean free path 
{L\ <C I) and the frequency is much greater then r , where r is the electron free 
flight time. As a result, the effect of this field on the charged particle is in a sense 
equivalent to the effect of a time-independent effective potential, that is a sequence 
of deep 'dynamic' quantum wells were the elelctrons are confined. The possible 
consequeces of this effect are also discussed and similarity with the classical Paul 
traps are emphasized. 

PACS: 73.20.D; 73.61; 85.30.V 
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1 Introduction 

To the author's knowledge, the term 'dynamic localisation' was coined by D.H.Dunlop 
and V.M.Kenkre in their well-known paper [1], where they considered charged particle 
motion on a liner chain of sites m (— oo < m < oo) under the combined action of 
a time-dependent electric field E(t) = Eosinut in the direction of the lattice and of 
nearest-neghbor coupling V. The Hamiltonian considered in [1] is of the form: 

oo oo 

H(t) = m >< m + 1 I + I m + 1 ><m |) — eE(t)d \ m >< m |, 
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where | m > represents Wannier state localised on the lattice site m, d is the lattice 
constatnt, V is the nearest-neigbbor intersite overlap integral and e is the particle charge. 
One of the main results of the paper [1] is that under certain conditions charged particle 
has to be localised in such a lattice, namely, the mean-square displasement < m > is 
to be bounded if the ratio of the field magnitude to the field frequency is a root of 
Jo, where J n , (n = 0, ±1, ±2, ...) is the ordinary Bessel functions in terms of which the 
correspoding exponentials are expanded. 

On the other hand, there is very similar but somewhat different type of charged 
particle localisation which also could be called 'dynamic'. It is the charged particle con- 
finement in a nonuniform rapidly oscillating field which for the first time was considered 
in the framework of classical physics by W. Paul and M. Raether [2]. This seminal work 
was the starting point for the development of something which was called 'Paul traps' 
and which later became one of the main tools for trapping a few and even single ions 
for the purposes of ions cooling and high-resolution spectroscopy. 

The point is that, as is known from calssical electrodynamics, there are no absolute 
maxima and minima of a potential in the electric field free of charges. As a consequence, 
the localisation of charged particles in this field is impossible if one means under locali- 
sation such a state where the particle with energy less than some definite value cannot 
leave the bound region under any initial conditions (Earenshow Theorem). However, as 
it was shown by W. Paul and M. Raether [2], the localisation of a particle in a non- 
uniform high-frequency electromagnetic field is possible. It turns out that under definite 
conditions the force acting on a particle can be represented by means of some effective 
time-independent potential. The force also does not depend on the sign of the particle 
charge and the particle can be localised in the space region where effective potential has 
its minimum. 

It is quite clear that consistent description of ions trapping and their subsequent 
cooling in Paul traps needs to be treated in the quantum mechanical framework. Such 
quantum-mechnical analysis of the particle motion in a rapidly oscillating field has been 
done by R. J. Cook, D.G. Shankland and A.L. Wells [3]. Having supposed the amplitude 
of the particle wave function slowly varying function of time (to compare with short 
time interval 2tt/uj, where u is the field frequency), and making average over time much 
greater than 2tt/u, the authors of [3] have shown that indeed, the particle motion in 
such rapidly oscillating field could be considered as if it occured in an effective time- 
independent potential V e ff = \/V (x) ■ \/V (x) / Am.uj 2 , where V(x) is the space-dependent 
part of the initial potential, m is the particle mass. This result is equivalent to the 

— » 

classical one, since the force is the negative gradient of the potential: / = — V V( x )- 

Somewhat later author also considered charged particle motion in a rapidly oscillating 
field in the quantum mechanical framework [4,5]. Contrary to [3], in [4,5] first, particle 
motion in a space-periodic high frequency field was considered and second, there was 
no any averaging over time. It turns out, that under certain conditions the solution of 
corresponding Scrodinger equation is asymtotically exact: if the spatial preiod is small 
enough and field frequency is sufficiently high, the effect of the field on the charged 
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particle is in a sense equivalent to the effect of the time-independent effective potential, 
that is a sequence of deep 'dynamic' potential wells where the electrons are confined. 
However, in [4,5] only one-dimensional model was considered, which is certainly not up 
to what could be required, since it is not relevant in full to the real world: in a real 
world charged particle always has a possibility to escape to other dimensions. Hence, 
if one would like to treat the results of [4,5] seriously, one should point out physically 
meaningful situation where such confinement could be observed. So, the aim of this 
paper is to consider more realistic model of such micro-scale dynamic localisation and 
to discuss its possible consequences. 

2 One-dimensional model 

Since the one-dimensional model plays an important role in further considerations, we 
discuss it here breafly putting the details of calculations into Appendix 1. 

So, let us have a charged particle (an electron, for definiteness) moving in a potential 
of the form V(x,t) = V(x) cos cut where V(x) = V(x + n£) = V coskox, and £ is the 
spatial period, n is the integer, /c = 2ti/£, Vq is the amplitude. Suppose the movement 
of electron is governed by the one-dimensional Schrodinger equation 



where e,m are the electron charge and mass, repectively. Suppose for a moment that 
only the potential-energy term were present on the right-hand side of (1). Then the 
solution of this equation would be 



This shows that dominant effect of the potential is to add an oscillating phase factor to 
the wave function tp. Hence, as it was suggested in [3], it is natural to look for a solution 
to Eq (1) of the form 




(1) 



ip(x, t) = ip(x, 0) exp(— iV(x) sm{u)t) /hw). 



ip(x,t) = ip(x, t) exp(— iV(x) sm(ut)/hu). 



(2) 



Substituting (2) and V(x) = V cosk x into (1) we have 



ih-£ = [H + H 1 (t) + H 2 {t)]<p, 



(3) 



where 



H = 




+ (4muj 2 ) ^^(eVokox) 2 sin 2 kox, 

eV k ) 2 . 2 , 
— sin k x cos loot, 



H 1 = 
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ih d 1 

H 2 = (ekoVosmkox — — h -eVok 2 . cosk x) smut. 

muj ox 2 

Now it is clear that averaging over time interval much greater than 2ti/uj, or in 

other words, over the time interval characterising changes of the amplitude tp(x,t) and 

substituting smut and cos2ut by their average values of 0, one would have the particle 

motion in effective potential 

= (eVofc ) 2 siii 2 k x = \/V(x) y V(x) 
eff 4mu> 2 Amu 2 

that is, the result of [3]. But we do not restrict the consideration only to that rough 
approximation; instead, we do carry out the subsequent analysis in two steps: first, we 
consider the equation with time- independent right-hand side (i.e., with the Hamiltonian 
Hq) and then, 'turning on' the terms Hi + H 2 , look at the consequences of it. We have 
then the equation 

dtp h 2 d 2 (eV k ) 2 

, % = ( -^^ + ^w- sm kox)(p - (4) 

Introducing the new variables x = k x and <p{x,t) = <p(x) exp(— ifi _1 e n f) we obtain 

-^ + {a n + 2q cos 2x)<p = 0. (5) 

Here a n = — /3(1 — e n /a), q = (5/2 = ma/tfk 2 , are dimensionless quantities; a = 
(eV k ) 2 /8muj 2 and E'fi £1X6 the quantities of the dimension of energy. 

Equation (5) is the Mathieu equation, the parameter a n , considered as a function 
of q, is the eigenvalue related to the corresponding Mathieu function. The Mathieu 
functions are the eigenfunctions of the Sturm-Liouville problem for the equation (5) and 
the boundary conditions either 

<£(0) = (p(-K) = 0, for se n (x, q) 

or 

T <*«> « 6 » 

Which of the two conditions and hence, which of the two Mathieu functions, se n (x, q) 
or ce n (x, q) is to be chosen, depends on the particular situation; it'll be discussed a little 
bit later. Now let us make some estimates. Suppose the characteristic length (spatial 
period) of the potential £ ~ 10" 4 cm, k = 2n x 10 4 cm _1 , cu ~ 10 10 Hz, V = IV, q = 
(l/8)(eV /huj) 2 = 2.89 • 10 9 . Hence, at so large values of q one can use the asymptote of 
the eigenvalues [6,7] 

a„ 2|g| + 2(2 + l)|g| 1/2 ; n = 0, 1, ... 

and for the energy values we have 

e n =(2n + l)\q\- 1 / 2 a. (7) 
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Note that the obtained 'spectrum' is equidistant with the distance between the levels 
As ntn -i = ehVoka / 4mw , which for the chosen values of paprameters is equal to Ae ni „_i = 
7.33' x 10- 15 erg = 4.5 meV . 

Let us consider now the solution of (3) with time-dependent Hi(t) + H 2 (t) term. In 
accordance with common quantum-mechanical rules, one should search for the solution 
to (3) in the form 

tp(x,t) =J2b n (t)<p n ex.p(-ite n /h) 

n 

Here (p n are the Mathieu functions from above and the time-dependent coefficients b n (t) 
obey the system of equations 

ihb n (t) = J2 H nm(t)b m (t), 

m 

where 

H nm (t) = HjU(t) + H®(t) = exp{iuj nm t) £ n (#i(x,t) + H 2 {x,t))(p m dx. 

Jo 

Here uj nm (£ n £m)- 

Depending on whether the functions are n or 27r-periodic, in the last integral the 
upper limit is n or In. 

Using the orthogonality and normalization of Mathieu functions, one can prove that 
Hnm = if and <p m are of different parity and H^l 7^ if <p n and y? m are of the same 
parity, that is both are odd or both are even. 

Just the opposite, if^ 7^ only if (p n and Cp m are of different pariry and hence, one 
can treat the terms Hj^ and Hj^ independently. It can be shown (see Appendix 1) 
that the system of equations which the coefficients b n are to obey, is asymtotically of 
the form 

ihb n (t) = £ H<U(t)b m (t) = £ M nm b m (t)exp(iuj nm t). (8) 

m m 

Here M nm is the real matrix explicit form of which is given in Appendix 1. Using an- 
other remarkable asymptotes of the coefficients of corresponding expansions representing 
Mathieu functions (for details see Appendix 1) and substituting C n = b n exp(i\ n ), where 
\ n are the parameters to be determined, one has 

%hC n = M nmC m exp[it(uj nm + A n - A m )]. (9) 

m 

As it is shown in Appendix 1, an appropriate choise of Aj removes the time dependence 
from if i(t), yeilding the solution C = C(0)e- ih ~ HM , where (7(0) = (Ci(0), C 2 (0), C N (0)) 
is the row matrix and iV is the maximum number of states (quantum levels) to be in- 
cluded in consideration. 

Since the matrix M nm is Hermitian, the solution of (9) finally has the form 

C(t) = diag[e^^*, eW ft > 7st , e^ n ^ nt ]C(0), 
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where diag[...] is the diagonal matrix, 7, are the eigenvalues of M, which are not nec- 
essary all to be different. So, despite the fact the Hamiltonian in (3) is explicitly time- 
dependent and the energy, generally speaking, is not conserved, the asymptotic proper- 
ties of Mathieu function caused by very large value of q and the small ratio of u/u ntn+ i 
cause the charged particles to be in the states with the energies e n and probabilities, 
corresponding to their initial distribution. 

The physical explanation of this consisits of the fact that, since u> <C w ni „ + i, the 
perturbation H(t) can be regarded as adiabatic. But as it is known from quantum 
mechanics, the adiabatic perturbation cannot cause the transitions between the states 
of the discrete spectrum [8]. The adiabaticity condition is that the changes of the 
interaction energy during a period of oscillations in quantum systems are much smaller 
than the absolute values of the energy differences between the corresponding states [8]: 

knm^ < n\H(t)\m > I < hu nm 

So, the existence of the if (t)-term in the Hamiltonian results only in the appearance of 
the new phase factor of the wave function stationary states and hence we can regard 
the charged particle as moving in some effective time-independent potential U e ff = 
2a sin 2 kox. Since a, having the dimension of energy, is much greater than the energy of 
the lowest eigenstates of this potential as well as of ksT, (ks is the Boltzmann constant) 
for any reasonable temperature, it is quite obvious that for the electrons this potential 
is the sequence of a deep 'dynamic' quantum wells where the electrons are confined or 
localised. However, it is also possible to provide more formal proof of electron localisation 
in such potential. The formal proof (see Appendix 2) is based on the observation that 
localisation also means the mean electron momentum in such a state should equal to 
zero. 



3 More realistic physical model: quantum wire 

Now let us proceed to the real world keeping however in mind and trying to preserve 
the benifits of one-dimensional model. How could we manage to do that ? Let us undo 
that 'knot' in a downright fachion, restricting the electron motion in other dimensions. 
Consider the quantum wire (see Fig.) containing two-dimensional electron gas (2DEG) 
under quantum Hall-effect conditions. It is known that under these conditions [9] the 
electron free path in 2DEG can be of the order of 100 nm and even somewhat greater. 
Suppose the total length L x of the quantum wire is L x ~ I — 100/im and suppose that 
on the surface of the wire thin insulator layer (its thickness will be estimated a little bit 
later) is deposited. If the quantum wire would be made of silicon, the insulator layer 
could be capping oxide. Suppose also that on the surface of the 'sandwich' there are 
the metal field electrodes in a periodic sequence with the spasing between the electrodes 
L\ ~ 5 x 10 _4 cm. To produce such grating is not a problem for the modern technology; 
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for example, in [10] the development of a microstructure with a distance between the 
field electrodes of lOOnm was reported. 

Now by means of suitably patterned gate electrodes, one could impose on a two- 
dimensional electron gas an artificial periodic potential of the form U(x,y,t) = (U\(x) ■ 
^(y)) cos out, where Ui(x + nLi) = U\{x). One could also suppose [^(y) to be weakly 
dependent of its argument, that is the first and second derivatives of [^(y) with respect 
to y(0 < y < L y ) to be very small. This assumption is quite clear and acceptable ; 
indeed, the potential applied to the gates is almost uniform along the metal strip that 
is, within the interval (0, L y ) and changes only at the very edges of the strip in the 
vicinities of y = and y = L y . So, one can conclude that its first and second derivatives 
are practically equal to zero almost evrywhere within the interval (0, L y ). As for U\(x), 
at the moment one cannot say anything escept the artificial potential imposed on 2DEG 
is a periodic one. But let us investigate more thoroughly the role of the insulator layer. 

Since the density of charges in the insulator is much smaller than in semiconductor, 
it is generally supposed to be zero and hence, one can consider the potential Q(x,y,z) 
within the dielectric as governed by the Laplace equation: 

dx 2 dy 2 dz 2 

If, as we supposed above, the electrodes are suitably patterned and the potential $ 
is weakly dependent on y, one can proceed from (10) to the equation 

<9 2 $ <9 2 $ _ 

dx 2 dz 2 

with the boundary condition 0,to) = Ui(x)cosut = CUi(x). Since $(x,0,to) is a 
periodic function with respect to x, &(x, z) could be expanded into a Fourier series of 
functions & n (x,z) = F n (z) cos{2nx/ L\), with the functions F n (z) obeying the equation 

d 2 F 



dz 2 



(4vrV/L 2 )F n , 



and hence, F n = Foexp(—z/z n ), where z n = Li/27m. Now it is obvious that the 
amplitude of the n-th harmonic decays exponentialy with the number n inreasing. The 
natural question is: how thick the insulator layer should be for one could neglect the 
higher space harmonics, restricting the consideration only to the first one ~ cosk x, 
where k = 2n/Li. If we suppose the amplitude of the second space harmonic is ten 
times smaller than the first one, we can calculate by means of the formulae above, 
that the insulator layer thickness should be about lOnm. This simple analysis shows 
quite clear that the insulator layer filters off the higher space harmonics of imposed 
potential and the electrons in the bulk of the semiconductor beneath the insulator are 
mainly affected by the lowest ones. The higher harmonics are strongly suppressed in the 
insulator and can be neglected if the insulator is sufficiently thick. 
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The next question to answer, is this: what the equation governs the electron motion 
in 2DEG under the circumstances considered above ? The answer may seem to be 
obvious, but in fact it requires a special care. 

It is well-known that the dynamics of electron in the semicnoductor conduction band 
under the time-independent external field can be described by an equation of the form 
[9]: 

[E e + { — ^ A)2 + U(r)Mv) = E9(r), 

where U(r) is the potential energy due to space-charge etc., A is the vector potential and 
m* is the effective mass. It should be noted that the wave function \l/(r) we calculate 
from this equation is not the true wave function but is the smoothed out version that 
does not show any rapid variations on the atomic scale (see for details [9]). 

Look again at the Fig. and suppose at first the constant potential is applied to the 
gate electrodes. Thus we suppose the electrones are free to propagate in x — y plane 
and are confined in z-direction. Usually at low temperatures with low carrier densities 
only the lowest subband corresponding to the confinement is occupied and the higher 
subbands do not play any significant role. We can then ignore the z-direction altogether 
and simply treat the semiconductor as a two dimensional system in x — y plane. 

By analogy with the last equation, suppose the dynamics of the electrons of lower 
subband in 2DEG under time- dependent external electric field applied to the gate elec- 
trodes, is governed by the Schrodinger-like equation of the form: 

Ti 2 

Here in accordance with the analysis presented above, U(x,y) = Ui(x)U 2 (y) where 
XJ\ = V cosk x , k = 2n/Li and U 2 (y) is supposed to be 



u 2 (y) 



Ci, < y < L y 
C 2 , y = 0,y = L y , 



where Ci( 2 ) are the constants (|Ci| > |C 2 |). Potenial U 2 (y) defined in this way means the 
jumps of it at the edges of the sample in y-direction, that is the potential discontinuity 
on the semiconductor-vacuum interface. The last one is equivalent to the Sommerfeld 
'rigid box' model of the solid state [11], for which boundary conditions for the electron 
wave funcrion x(y) are: x(0) = x(L y ) = 0. 

Searching for the solution to (11) in the form \l/(r) = x(y)(/9(a;, t) = exp(ik y y)ip(x,t), 
one gets the variable separated and together with the boundary conditions, the next sys- 
tem of eqautions holds: 

h 2 d 2 

2 x(y) = z y x(y) 



2m* dy 1 

., dip r h 2 d 2 TT . . . 
lh m= [ 2m-*¥x- + Ul{x)C ° SUjt]iP 
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Here 

£y = 2rt k »' ky = Tj 1 ' n = °' ±1 ' ±2 "' 

With C/i(x) = eV^ cos fco^, the second equation becomes the same as that has been 
considered in previous section. It is clear that in order the approach discussed above to 
be valid, in particular, for the transformation (2) can be used, the frequency applied to 
the gate electrodes should be greater than the inverse time which takes for the electrones 
to fly between two successive gates. The last one can be estimated as follows. Let it 
be the electron concentration in 2DEG equal n s = 5 • I0 u cm~ 2 and the Fermi velocity 
is of the order ^~3x 10 7 cm • s -1 . Then the flight time r = L\jvp pa (5/3) • 10 _11 s 
and hence, r" 1 w (3/5) • 10 11 Hz. It means the condition u ^> r 1 could be fulfilled if u 
would be at least of 10 12 Hz. However, it is worthy to note that in accordance with the 
formula for A£„ ;n _ 1 , increasing the distance between the gates by a factor 2, one could 
reduce the friequency of the field applied to the gates by factor 4, retaining the spacing 
between the levels intact. Thus, if L\ ~ 10~ 3 cm, the frequency can be of ~ 2.5 • 10 11 .Hz. 

Also we should discuss which of the two boundary conditions (6) for the Mathieu 
equation has to be chosen. Having in mind possible experimental checking up the pro- 
posed model, it seems natural to consider quantum wire connected to the current leads. 
Then the second of the two conditions (6) is to be more appropriate, since the first one 
means there is no any current at the quantum wire-lead interfaces. 

To complete analysis, make again some estimates. Suppose L\ = 10 _3 cm, oo ~ 
2.5 • 10 n Hz,m* pa 0.1 • m e (m e is the free electron mass), Vq = 10V. Then, q pa 4.628 • 
10 8 , Ae niri _! pa 4.5meV and Ae n>n -i/hu! ~ 28 . As a result, we have come to the next 
conclusion: if in a quantum wire with the periodic gate electrodes considered above and 
under quantum Hall effect conditions, the spacing is small enough and the frequency of 
the field applied to the gates is sufficiently high, the effect of the field on the charged 
paricles is in a sense equivalent to the effect of a time-independent effective potential. 
The 'potential relief is a sequence of deep 'dynamic' quantum wells of the depth 2a and 
the spectrum defined by the formula (7). The physical explanation of this at first sight 
paradoxical, fact is that, since Hcu <C Ae ni „_i, the perturbation H(t) can be regarded as 
adiabatic and hence, it cannot cause the transitions between e n - states. 

Also it is curious that, despite of the initial potential U(x) = Vocoskox does not 
look like the harmonic one of the standard quantum mechanical textbooks, the resultant 
specrtum turns out to be equidistant just like for harmonic potential ~ x 2 it is. This fact 
also can be easily explained. Indeed, at chosen values of the structure parameters, the 
main dimensionless parameter q which defines all the physics, is of the order of 4.628- 10 8 
and hence, the effective time-independent potential U e ff = 2asin 2 koX for small values 
of x looks like the harmonic one. 
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4 Discussion and conclusion 



In the paper the dynamic localisation of charged particle (electron) in a quantum wire 
under the external non-uniform time-dependent field is considered. The electrons are 
trapped in a deep 'dynamic' quantum wells (stricktly speaking, quantum dots, since 
the electron 'energy spectrum' turns out to be completely discrete in there) which are 
the result of specific features of the potential imposed on 2D electron gas: the scale of 
spatial nonuniformity is much smaller then the electron mean free path (L x <C I) and 
the frequency is much greater then r _1 , where r is the electron free flight time. As a 
result, the effect of this field on the charged particle is in a sense equivalent to the effect 
of a time-independent effective potential, that is a sequence of deep 'dynamic' quantum 
wells were the elelctrons are confined. 

Certainly the next question immediately arises. Whatever large, the dimensionless 
parameter q is not however, infinitely large. But mathematically it is just not we need: 
we should have q — > oo, in order to neglect some part of the perturbation. So, what are 
the consequencies of the q is finite ? 

It is quite clear that the picture sketched above is just an approximation. Since q 
though large, is finite, neglected part of the perturbation would cause the transitions 
between the levels. However, the rate of the transitions is small compared to the inverse 
field frequency. There are should be about 20 or 30 of field oscillations for the transition 
probability reaches the unit. The complete analysis of charge carriers transport in the 
'dynamic' quantum wells under such transitions is beyond the scope of the paper; we plan 
to discuss it in the next publication. Here we note only that the localisation considered 
above seems to lead on one hand, to the current suppression in a quantum wire and 
on the other, to the emerging of current peaks on its current-voltage characteristic (so 
called, I-V curve). These peaks are due to electron delocalisation, or escaping from 
the wells by means of resonant tunnelling which seems to be quite possible, since the 
tunnelling time (do not confuse it with the tunnelling rate which is proportional to the 
inverse tunnelling probability) is much smaller then the rate of transtions between the 
e n -states. Indeed, omitting the subtle question of prooper definition of tunnelling time 
(see discussion on that subject in [12]), we could estimate it roughly as r« ~ Sb r /v g , 
where Sb r is the width of the 'barrier' between two successive dynamic wells and v g - 
group velocity which could be supposed equal to Fermi velocity. Supposing s& r to be 
approximately ten times smaller than L\ (which seems quite reasonable), we get then 
Tu 200 or even 300 times smaller then the transmission rates. Hence we conclude that 
for tunnelling electrons the 'potential relief they encounter on their way indeed looks 
like the static one. 

Putting it into plain words, one can say that since the electron life time in e n - states 
is long enough (the transition rates as we remember, slow down), no wonder that just 
for these states probability of tunnelling increses. 

Suppose we have a sequence of barriers on which partciles are incident and suppose 
the energy of particles is smaller than the height of barriers. Then the barriers are 



10 



practically impenetrable for almost all particle energies; however, for certain discrete 
energies and respective energy level widths particle can pass through the barriers without 
any reflection. This is the resonant tunnelling phenomenon [13,14]. It is quite remarkable 
that from the mathematical point of view the transparency of the barriers for these 
discrete energies does not depend on the width of the barriers if the barriers are identical 
[13]. 

Return now to the stationary equation (4) and rewrite it in the form: 

^J" + (e - U eff (x))^ (12) 

where 

U eff (x) = 



2a sin 2 kox, x e [0, L x ] 
0, x < 0, x > L x . 



The tunnelling problem consists of finding all the solutions of (12) parametrically 
dependent on e G (0, oo) and behaving as 



1> 



exp(ipx) + r(p) exp(-ipx) x < 0, 
t(p) exp(ipx) x > L x ; 



(the case of a normalized particle beam incident from the left); p 2 = e; t(p), r{p) stand 
for transition and refelection amplitudes, respectively. 

We introduce now the outgoing propagator G + (x,x';p) and define the electron wave 
function along the internal region < x < L x by ip(p,x) = 2ipG + (0 , x; p) . Then 
the transmission amplitude is t{p) = 2ipG + (0, L x ; p) exp(— iL x ). Near an isolated pole 
Pn = $n — iln it is possible to write the propagator in the form [13]: 

Here (p n (x) are the eigenfunctions of the Hamiltonian in the left-hand side of (12), 
p\ = e' n — iT n /2,e' n = 5 2 — 7 2 ,T n = 4(5„7„. Using (13) and taking into account that 
5 n In, for the transmission coefnceints |t(p)| 2 we have 

\t(p)\ 2 p2 ^ n ^ 2 ^ Lx ^ 2 



Let us notice that in accordance with [13], 7„ = (|(/? n (0)| 2 + \ip n (L x )\ 2 )/2I, where 
I m 1. It is also clear from the formula for |t(p)| 2 why we choose the second of the two 
boundary conditions (6). 

If the potential in (12) possesses only two maxima (consists of two barriers), i.e. 
L x = 2ii/k , then, since the barriers are symmetric, i.e. |<^ n (0)| 2 = \(p n (L x )\ 2 , we have 

lim 
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Now, let the number of identical barriers be m > 2; then as it was proved in [14], if 
any arbitrary pair of barriers is completely transparent at some energy e = e n , the whole 
structure consisted of m > 2 replicas is also completely transparent at this energy. 

Thus it may be said that just like in the static quantum wells, in the sequence of 
'dynamic' ones, the electron resonant tunnelling also seems possible because the barrier 
transmition coefficient has maxima if the electron energy e — > e n . 

To our mind, such 'dynamic' quantum wells would have even some advantages in 
comparison with static semiconductor heterostructures, since their 'spectrum' could be 
controlled by the external high-frequency electric field applied to the gate electrodes. 

In previous section we estimated the field frequecny uo as to be equal ~ 2.5 x lQ ll Hz. 
This value is determined mainly by the electron mean free path which under quantum 
Hall effect conditions can be about 100/im and hardly could be greater. To apply the 
field of such frequecny to the microstructure, whatever possible, is difficult problem to 
solve and it is interesting to know whether it is possible to reduce the frequency. We 
plan to consider this really intriguing problem in the next publication, here instead, we 
discuss another interesting question: how the effect considerd in the paper relates to the 
quantization of charge carrier energy spectrum in the uniform electric field, that is to 
Wannier-Stark effect. 

In 1960, G.H. Wannier [15] studied electronic states in the presence of a uniform 
electric field and found that eigenstates are localised along the direction of the electric 
field and have quantized energy levels e n = neEd, where n is an integer, E the electric 
field and d the lattice period along the electric field. 

The charge carriers in a crystal in the external electric field, in order to pass from 
the domain of the lower potential to that one of the higher potential, should gain an 
additional energy. Hence, the electrons in a crystal in an external electric field are to 
be localised in a domain with the characteristic length of the order of 2A s /eE, where 
A s is the subband width. As a result, the electron energy spectrum in a pure crystal 
whithout impurities would resemble the ladder which is called Stark, or Wannier-Stark 
ladder. Therefore, the electron dynamics along the field direction, say x-direction, can 
be characterised by the ^-momentum only for those values of E, for which localisation 
length are much greater than d. This field range, 2A S ^> eEd is called classical. If the 
field increases, the situation may come into being when the localisation length becomes 
of the order of d. It is clear that the energy dependence on the qausi-momemtum k x 
is no longer parabolic, the electron turns out to be trapped in the potential well and 
its dynamics practically does not depend any further on the crystal lattice parameters. 
The corresponding field range, 2A S > eEd is called quantum. 

The localised Stark ladder states are associated with Bloch oscillations [16]. The 
frequency of Bloch oscillations is ujb = eEd/h. In order to observe the Bloch oscillations, 
one should have the next condition fulfilled: cubt p > 1, where r p is the momentum 
relaxation time. Since usually in the 3D semiconductors r p is of the order of 10~ 13 s, the 
field should be about 2 • lQ b Vcm~ l in order to satisfy this condition. Usually even the 
lower fields cause the electric breakdown and that is why the Bloch oscillations were not 
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observed in bulk semicnoductors. 

There have been also a number of other controversial arguments on the Stark ladder 
states. Some claimed [17] that such ladder states cannot exist because they will decay 
rapidly due to Zener tunnelling caused by interband mixing. Some others claimed that 
such interband effects are totally absent and no Zener tunnelling is observable [18]. 
Others claimed that Zener tunnelling is possible in Stark ladders [19,20]. 

A number of experiments has been done in order to observe a Stark ladder in bulk ma- 
terials, but no clear evidense was found. Only weak oscillations of conductivity observed 
in ZnS was attributed to the hopping motion of electrons between Stark ladder states 
[21]. The failure of observation is attributed to scattering by impurities or phonons, 
which prevent the accelaration of electrons up to a Brillouin-zone edge. However, in 
superlattices electrons can be accelerated easily to the edge of the Brillouin zone before 
they are scattered because of the large lattice constant. Since then, a number of ex- 
periments on Wannier-Stark effects have been performed using various techniques (see 
[22]). 

It is clear that the effect discussed in the paper though different, resembles to some ex- 
tent Wannier-Stark localisation , because it is the electron localisation in high-frequency 
non-uniform electric field. Since the localisation length is much greater than d and 
the momentum relaxation time in quantum wires under quantum Hall effect conditions 
is much greater than in bulk semiconductor, one can hope that these states could be 
observable just like Wannier-Stark states are observable in the superlattices. 

For example, it is known that magnetophonon resonance (MPR) is observable in 
the superlattices [23]. The MPR is due to the oscillating behaviour of the electron 
density of states arising as a result of the Landau quatization. The magnetophonon 
resonance appears every time when the phonon frequency uulo is equal to the cyclotron 
frequency uj c in a magnetic field multiplied by small integer n, ujlo — nu c ,n = 1,2, .. 
where u c = eB/m*. Thus, it seems probable, that in case of the localisation effect 
discussed here, the current peaks on the I — V characteristics of the structure also 
could be observed when the optical phonon frequency uj lo would be equal to Ae nm /h. 
Here e n -states would play the same role as the Landau levels do in MPR and the last 
phenomenon, if observed, could be called electrophonon resonance in the high-frequency 
non-uniform electric field. 

As the final remark, it is worthwhile to return again to the beginning of the paper in 
order to compare the effect discussed here with the Paul trapping effect and emphasize 
their striking similarity. 

Indeed, the ion motion in radio-frequency non-uniform electric field to the first ap- 
proximation can be represented as the sum of comparatively slow motion r(t) in some 
effective potential and the rapid oscillations with small amplitude rj(t) near the local 
equilibrium r (t): 

r(t) = r (t) + r}(r ) cos S2 t, 
where VLq is the applied field frequency and the amplitude of small oscillations r)(t) is 
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determined by the trap electric field: 



r](r ) = eE{v )/MQl 
here M is the ion mass and the effective potential is of the form [24,25]: 




Parameter A is defined as A = U /(pl + 2zq), U is the amplitude of the field and p , zo 
are the cylindrical coordinates of r . In this effective potential an ion oscillates along 
z- axis with the frequency uj and in the x — y plane with the frequency uj p = a>/2, where 

cu = 2V2eA/MQ . 

More thorough analysis [25] has shown that the oscillations in the effective potential 
<& e ff are stable if the parameter is sufficiently large. This analysis also showed the 
ion oscillation spectrum consists of the infinite set of discrete frequencies nfl ±u,n = 
0,±1,±2.... 

It is easily seen that our condition uj ^> r -1 corresonds to the stability parameter 
introduced in [25] and the states e n correspond to the spectrum of ion oscillations in the 
Paul trap. 
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6 Appendix 1 

Here in Appendix 1 we derive the system of equations (8) to which the ^-coefficients 
are obeyed. Let the Mathieu functions </?° and (p Q m be of the form (see [6]) 

oo 

$n = ce 2n (x,-q) = (-l) n ^(-l) r 4 2 r n) cos2r£, 

r=0 

oo 

V m = ce 2n (x,-q) = (-l)«]T(-l)VL£ m ) C os2r£. 

r=0 

Using the orthogonality of the functions one can prove that = af nm cos2cut, where 

1 oo oo 
^ r =0 s=o 

= --{A {2n) A^ m) + Af n) A { 2 2m) + ... + A [ Q n) A { 2 2m) + A { 2 2n) A {2m) + ....). 

The evaluation of the sum in the brackets requires special care, but one can note 
that the asymptotic behaviour of A^ at q — > oo does not depend on the superscript. 

For example (see [6]): lim^^Mc^) = -2, lim^^M^) = (~l) r 2 at 
arbitrary n. 

Using the recurrence formula for the function ctin 

a 2n A {2n) - qA {2n) = 0, 

one can demonstrate that for q = 1600 the ratio A 2 /Aq = —1.95 for the function ceo 
and -1.75 for the function ce 2 . At the q = 4.628 • 10 s we have A 2 /A = -1.999907 for 
ce and —1.99972 for ce 2 , respectively. Hence we conclude that 

fnm ~ - ^^2^^-2(^-1) 

s 

and 

H nrL = a fnm cos 2ut exp(iu nm t) = ^af nm {exp i(u nm + 2u)t + exp i(uj nm - 2w)t). 
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Note that at the chosen values of parameters 2cj <C uj nm , and hence 

w M nm exp(iu nm t), M nm = af nm /2. 

Suppose now that 

<p n = ce 2n+1 (x, -q) = (-l) n E?lo(-l) r s£S 1) cos(2r + l)x, 
Vm = ce 2m (x, -q). 

Then H^) n = fig nm sin out, where \i = —i2hekl(yo/mu) and 

oo 

-X^ r A ( 2m ) I R (2n+1) , R (2n+1) x 
i/nm — 2^ ' ^2r l-°2r+l -°2r+3 )• 

r=l 

The last series is convergent, since coefficients A^™\B^V^ are of the order of r~ 2 
at r — > oo. It is noteworthy that the coefficients -B^r+i^ a l so do n °t depend on the 
superscript at q — > oo: lim^oo B^"^ 1 V-^i 2n+ ^ = ( — l) r - So, if^ has the form 

H nl = ^^(expi^nm + - expi(cj nm - w)t) 

and for g — > oo and uo/uo n ^ n+ i <C 1, i/^ — > . Finally we have 

exp(iuj nm t)b m (t) , (15) 

that is, Eq (8). Owing to the absolute convergence of the series (14), one can conclude 
that f nm does not exceed some finite quantity. Also it is known [6] that A^^ -coefficients 
possess another remarkable asymptote, namely A^^ — ■> for n — > oo and n > q. Hence, 
we can restrict the number of states to consider by a finite number N. Substituting 
C n = b n exp(i\ n ) where \ n are the parameters to be determined, one gets 

ihC n = M nm exp it(uj nm + A n - \ m )C m . 

m 

Thus, any choice of Aj such that for all n, m 

Unm + A n - A m = (16) 

removes the time dependence from Hi(t), yelding the solution C — (7(0) exp(— ih^Mt), 
where C = (Ci(0), C 2 (0), C N (0)) is the row matrix consisting of coefficients at the 
initial time t — 0. 

It is obvious that in general (16) includes up to (l/2)iV(JV — 1) equations in the AT 
unknowns Aj to be solved simultaneously, so that the consistent choice of the set {Aj} is 
not always possible. However, we have to recollect that our 'spectrum' e n is equidistant 
and in case of iV states the set {uj nm } consists of N — 1 unequal terms only. So, if we 
believe any of the Aj is equal to some definite value (zero, for example) we obtain N — 1 
simultaneous equations in N — 1 unknowns. Thus, as it is stated in Sec.l, we have 

C{t) = diagfe^ 71 *, e^ R > 7 «*, e {l/nhnt ]C{0), 

where diag[...] is the diagonal matrix and 7 n are the eigenvalues of M- matrix . 
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7 Appendix 2 

Here we prove the mean electron momentum is equal to zero in the potential U e ff = 
2a sin 2 kox, using the the Wigner function method [26]. 

Soppose the small subsystem described by the Hamiltonian 

d2 o • 2, 
H = -— — -— + 2a sin fc a: 

2m* dx z 

is in the state of thermal equilibrium with the environment characterised by the tem- 
perature T. Thus, the density matrix of the subsystem is of the form [8] : 

p(x,x') = Z _1 ^^*(x)^ n (x') exp(-e n /k b T), 

n 

where statistical sum 

Z = J2eM-tn/k B T) = Sp[(-H/k B T)} 

n 

allows the density matrix to obey the normalization condition Sp[p] = 1. 

The diagonal element p(x, x) = P(x) is a probabilty for a particle to be at the point 
with the coordinate x, while p(p,p) = Pip) is a probablility for a particle to have the 
moment p. 

In the framework of classical statistical mechanics there can be introduced probability 
density function with the properties: 



p ip) = J f(p,x)dx, P(x) = J f(p,x) 



dp 
2ttK' 



It turns out that in quantum mechanics the role of probability density function plays 
Wigner function fwip, x ) defined as [26]: 

fw(p, x) = J p(x + r]/2,x- rj/2) exp((-i/h)pr])dr]. (17) 

Then, for a function hip, x) which is the function only of p, or only of x, the next relation 
holds: 

/dp 
f w (p,x)h(p,x) 1 ^dx, 

where < .. > stands for mean value. Thus, in order to calculate < p >, we need 
in accordance with (17), to have the expression for the density matrix p(x,x). With 
Hamiltonian considered, the equation for density matrix reads 

ih^pix, x') = {H- H')p(x, x'; t). (18) 
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Since the Hamiltonian is time-independent, one can search for the stationary solution 
to this equation. It is easily seen, that if one searches for the solution in the form 

p(x, x') = C exp(—f(x + x', x — x')), 

it is possible, by the proper choise of the auxiliary function f(x + Ob . Ob Qlj ) to subtract 
the potential U e ff, sutisfying in this way the equation (18). 
One can check it up that the proper choise is this: 

171* '' 

p(x, x 1 ) = C exp[ ( — - ) 2 sin(A;o(a; — x')) — 2aa cos(k (x + x'))], (19) 

a 2n.ko 

where a is the arbitrary constant and C is the normalizing constant defined by the usual 
condition 



/ p(x, x')dx = 1. 
Jo 



Supposing a = ip, upon expanding the exponential of the trigonometric function in 
terms of Bessel functions [27] : 

oo 

exp(iysmtf) = £ J n (y)e in& , 

n=—oo 

where y = (m* / p)(l/2hk ) 2 , d = k r] and using (17), (19), by means of direct calculation 
one arrives at < p >= 0. 
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